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Abstract. This work evaluates the model dependence of the electric and Coulomb quadrupole amplitudes (E2, C2) in the 
predominantly Ml (magnetic dipole-quark spin flip) y*N — » A transition. Both the model-to-model dependence and the 
intrinsic model uncertainties are evaluated and found to be comparable to each other and no larger than the experimental 
errors. It is confirmed that the quadrupole amplitudes have been accurately measured indicating significant non-zero angular 
momentum components in the proton and A. 

PACS: 13.60.Le, 13.40.Gp, 14.20.Gk 



PHYSICS MOTIVATION 

Experimental confirmation of the presence of non-spherical hadron amplitudes (i.e. d states in quark models or p 
wave 7Z-N states) is fundamental and has been the subject of intense experimental and theoretical interest (for reviews 
see 01 |H [H). This effort has focused on the measurement of the electric and Coulomb quadrupole amplitudes (E2, 
C2) in the predominantly Ml (magnetic dipole-quark spin flip) y*N — > A transition. Since the proton has spin 1/2, 
no quadrupole moment can be measured. However, the A has spin 3/2 so the Y*N — * A reaction can be studied 
for quadrupole amplitudes in the nucleon and A. Due to spin and parity conservation in the y*N(J K = l/2 + ) — > 
A(J n = 3/2 + ) reaction, only three multipoles can contribute to the transition: the magnetic dipole (Ml), the electric 
quadrupole (E2), and the Coulomb quadrupole (C2) photon absorption multipoles. The corresponding resonant pion 
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production multipoles are M/ , E l+ , and S l+ . The relative quadrupole to dipole ratios are EMR=Re(£' 1 ^ /Mj^_ ) 

3/2 3/2 

and CMR=Re(5j^ /M/ ). In the quark model, the non-spherical amplitudes in the nucleon and A are caused by the 
non-central, tensor interaction between quarks Jj]. However, the magnitudes of this effect for the predicted E2 and C2 
amplitudes[5] are at least an order of magnitude too small to explain the experimental results and even the dominant 
Ml matrix element is ~ 30% low lHH. A likely cause of these dynamical shortcomings is that the quark model does 
not respect chiral symmetry, whose spontaneous breaking leads to strong emission of virtual pions (Nambu-Goldstone 
Bosons)[ 3]. These couple to nucleons as rj ■ p where a is the nucleon spin, and p is the pion momentum. The coupling 
is strong in the p wave and mixes in non-zero angular momentum components. 

However, the multipoles are not observables and must be extracted from the measured cross sections. The five-fold 
differential cross section for the p(e,e'p)n° reaction is written as five two-fold differential cross sections with an 
explicit 0* dependence as |0] 
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= r(or + £Gl + vltGlt cos <f)* + £Ott cos 20* + hp e v LT i<5 LT t sin0*) (1) 



dSlfdEfdO. 



where e is the transverse polarization of the virtual photon, vlt = \/2e(l + e), v LT i = y/2e{\ — e), T is the virtual 
photon flux, 0* is the pion center of mass azimuthal angle with respect to the electron scattering plane, h is the 
electron helicity, and p e is the magnitude of the electron longitudinal polarization. The virtual photon differential 
cross sections (07-, Ol, (Jlt, GtTiGlt 1 ) are ai l functions of the center of mass energy W, the four momentum transfer 
squared Q 2 , and the pion center of mass polar angle 0*^ (measured from the momentum transfer direction). They are 
bilinear combinations of the multipoles 16[] . 



RESONANT MULTIPOLE FITTING 



The current experiments do not have sufficient polarization data to perform a 

model independent multipole analysis and must rely upon models for the non-resonant (background) amplitudes. The 
standard procedure to extract the multipoles is to use the models to fit the data. Their background terms are unaltered 
and the three isospin =3/2 resonance multipoles 

R^ = Mli 2 ,E% 2 ,S% 2 (2) 

are fit to the data. Specifically we introduced multiplicative factors, A (/?,-) for the 7 = 3/2 multipoles so that the phase, 
and hence unitarity, is preserved. For the charge channels with a proton target and an outgoing neutral pion (e.g. 
7+ p -> 7t Q p): 

R l {K°p)=R)l 2 + 2 -l{R i )R]' 2 (3) 

where Ri represents any of the three photo-pion multipoles Mi + ,E\ + ,Si + for the final charge state or in the isospin 
1/2 and 3/2 states. 

Three parameter, resonant multipole fits were performed on data taken at Q 2 = 0.060 1 18] and Q 2 = 0.126 (GeV/c) 2 
l9 [ l l0LlllLll2Lll9ll. one Q 2 value at a time, using four representative calculations: the phenomenological MAID 2003 
1 20h and SAID 12 lll models. and the dynamical models of Sato-Lee i22ll and DMT 1I23I1 . The fits are presented in terms 
of Mj{ 2 , EMR = E2/M1 = Re(#i+/Mi + ), and CMR = C2/M1 = Re(5i+/Mi+). At least one multipole is expressed in 
absolute terms rather than as a ratio because some models can give accurate predictions for ratios but not for absolute 
sizes. Figure Q] shows our new, low Q 2 data along with several model predictions before and after the three resonant 
parameter fitting. The convergence is rather significant. Only the data points in the top three plots were included in 
the fits and yet the parallel cross section as a function of the center of mass energy W converged nicely. Note that 
since the Sato-Lee model does not include higher resonances, it was not expected to fit the data well at higher W 
explaining the deviation observed in Fig.Q] Also, as expected, the a LT i curves did not converge since this time reversal 
odd observable [24] is primarily sensitive to background amplitudes and the fit is only for resonant amplitudes. 

The bottom half of FigureQ]shows the "spherical" calculated curves when the resonant quadrupole amplitudes {E^ 2 
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in <7o an d S l+ in Gtl, see Appendix for multipole expansions of the observables) are set equal to zero. The difference 
between the spherical and full curves shows the sensitivity of these cross sections to the quadrupole amplitudes and 
demonstrates the basis of the measurement of the E\ + and S\+ multipoles. The small spread in the spherical curves 
indicates their sensitivity to the model dependence of the background amplitudes. Figure [2] shows the values for the 
EMR and CMR for the four models before and after fitting. It is seen that there is a very strong convergence of these 
values after fitting. We have quoted the average value of these parameters as the measured value and are using the 
RMS deviation to estimate the model-to-model error lll2l Il8ll since these four models are sufficiently different to have 
a reasonable estimate of the present state of model dependence of the multipoles. At the present time the model-to- 
model and experimental errors are approximately equal. 

In a way what we are observing is the fact that the electro-pion production process shows us two separate faces, 
depending on the observable and on the center of mass energy W that we choose. The best way to extract the three 
resonant amplitudes is to measure the time reversal even observables (ob, Orr, Glt) EH] at or near the resonance 
energy W — 1232 MeV. On the other hand, the best way to test the model calculations is to examine time reversal odd 
observables such as a LT i [24] right on resonance. In addition, we also have off resonance data. These are sensitive to 
both the shape and phase of the M\ + multipole and also the background amplitudes. The yp — > n + n charge channel 
is also more sensitive to the background amplitudes particularly the 1=1/2 amplitudes. Such background sensitive 
data in combination with model studies are essential if the field is to progress to the stage where the model errors are 
significantly smaller than the experimental ones. 





FIGURE 1. Q 2 = 0.060 (GeV/c) 2 data with model predictions before fitting (top panels) and after the three resonant parameter 
fit (bottom panels) along with the EFT predictions from Pascalutsa and Venderhaeghen (PV) (2^1 . Note the convergence of the 
models except for the background sensitive au' points. Data from 1 18] and include the statistical and systematic errors. The lines 
with dots on them are the fitted models with the E\ + and S\+ quadrupole terms set to zero. The models are MAID 2003 [20], SAID 
(H, DMT |H and Sato-Lee (SL) (H. 
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FIGURE 2. Q 2 = 0.060 (GeV/c) 2 extracted EMR and CMR before and after fitting. The light error band is the average of the 
fitting errors and the darker band is the RMS deviation of the models added in quadrature. The models are MAID 2003 1 20], SAID 
l2lll . DMT 123] and Sato-Lee t22ll and the effective field theory calculations of Pascalutsa and Vanderhaeghen (PV) l25ll and Gail 
and Hemmert (GH) (H. The convergence of the models after fitting is clear. 



INTRINSIC MODEL ERRORS IN DETERMINATION OF THE RESONANT 

MULTIPOLES 

Beyond Three Parameter Fits: Including Background Multipoles 

This work expands the three resonant parameter fits to include the influence of the background multipoles on 
the resonant amplitudes derived from fitting the experimental data. In this way we will be able to make reasonable 
estimates of the intrinsic model errors due to uncertainties in the background multipoles and to see if this leads to any 
suggestions to reduce them. First, we include the remaining s and p wave multipoles: Eq + ,Lq + , M\-,L\-. Next, we 
estimate the influence of the higher partial waves using the CGLN invariant amplitudes F(, i = 1, . . . ,6 Og]. We introduce 
a new combination of higher order multipoles we call Fj. These combinations show that many small multipoles can 
have a cumulative effect as will be seen shortly. The fjs are varied using a scaling factor A, as 

F t = F, s&p + X,Fi 

Fi = £ , o + + 3Mi + cos0 + 3£'i + cos0 + AiJY 

F 2 = Mi_+2Mi+ + A 2 ?2 

F 3 = 3(E l+ -M 1+ )+X 3 T 3 

F4 = X4F4 

F5 = Lq + +6Li + cos 6 + A5F5 

F 6 = Li_ — 2Li + + X(,F( > . 

Next, we allow the 1—1/2 part of the charge channel multipole to vary in a way similar to the / 

R i (n Q p)=X(Rf)Ry 2 + X(R l )^ /2 . 



(4) 

= 3/2 part. 

(5) 




This new fitting procedure introduces thirteen background amplitudes, too many to determine with the available 
data. So, we set out to determine how much the resonant parameters are affected by the uncertainties in the background 
amplitudes. We try to quantify this effect and to determine which parameters have a strong effect and which are 
correlated. 

First, we look for parameters which are correlated with the resonant parameters. Figure [3] shows some examples of 
negative, positive and no correlation. Fitting parameters are plotted on each axis and the ellipse indicates the region 
where 68% (1 a) of the fits are expected to fall if many similar data sets are fit. The ellipse with axes close to the x 
and y axes shows no correlation. However, the other ellipses are rotated indicating that as one parameter tends one 
direction, the other parameter tends to go with it or away from it. This indicates a correlation between the parameters. 
While error ellipse plots are useful in a qualitative way, they are difficult to use in a quantitative manner. Changes in 
the scale of the axes will change the angle of the ellipse hiding or exaggerating correlations. Therefore, we use the 
correlation coefficient, r, to indicate the level of correlation between two parameters: 
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are used 02711 . r varies from -1 to 1 and is insensitive to the parameter scales since the scale factor for each parameter 
cancels in the ratio. Figure|4]shows the various correlation coefficients for each background amplitude with E^,m\'^, 
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and Ly + using MAID 2003 and combined Bates and Mainz data at Q = 0.126 (GeV/c) . The square of r indicates 
how much of the variance is explained by a linear relation between the two variables. A rule of thumb is that two 
variables are correlated if r 2 > 0.5 and uncorrected if r 2 < 0.1. This then leads to the ranges in r. 0.7 < |r| < 1.0 = 
large correlation, 0.3 < \r\ < 0.7 = medium correlation, 0.0 < |r| < 0.3 = small correlation. 

The next check is for sensitivity. If a parameter is large, zeroing it out should affect the % 2 by a large amount. For 
example, when Lq+ is turned off, the model predictions change noticeably (see Fig. 0. Other background terms can 
have significant effects as well like the FjS and the remaining s and p do in Fig. [6] Table Q] shows the ^ 2 /d.o.f. that 
results from turning off the various background amplitudes in the MAID 2003 model. It also indicates how strongly 
the amplitude was correlated with any of the resonant multipoles. 

In Figs. |5] and |6l a combined cross sectionlB IToL [l2ll Oei = Oo(9* q ) + 0rr(0;L) — OQ(6* q — n) is shown. In this 
linear combination the dominant M\ + multipole contribution cancels out and shows the effect of the smaller isi-i- 
quadrupole contribution. (See Appendix for the expansion of the observables in terms of multipoles.) 

Using the criteria of correlation (\r\ > 0.7) and sensitivity (x 2 increases by 50% upon removal), several parameters 
were identified as significant using the Q 2 = 0.126 (GeV/c) 2 data set. Those multipoles are shown in Tableland 
include two of the s and p multipoles, two of the isospin = 1/2 multipoles and three Fj terms. Looking at Fig. [7] all 
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FIGURE 4. Correlation coefficients for background amplitudes with the resonant amplitudes E l+ ,M l+ ,L l+ using MAID 2003 

for Mainz and Bates Q 2 = 0.126 (GeV/c) 2 data. The central region of each plot is for small correlation with the next region being 
medium correlation and the next being large correlation. Note: Flb=Fi, F2b=F2, etc. 



TABLE 1. Gives £ 2 /d.o.f. resulting 
from turning off the corresponding 
background parameter in MAID 2003 
arranged from most to least sensitive us- 
ing the combined Bates and Mainz Q 2 = 
0.126 (GeV/c) 2 data. The type face in- 
dicates the level of correlation with any 
of the three resonant multipoles: Large 
Correlation Medium or Small Correla- 
tion. 

Extra Par. ^ 2 /d.o.f 

7.42 

6.09 
5.60 

5.59 
4.10 

3.31 
1.85 

1.85 
1.80 

1.72 
1.55 
1.24 
1.21 
1.17 
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FIGURE 5. Sensitivity to L + in MAID 2003 using data at Q 2 = 0. 126 (GeV/c) 2 . The solid curve is the full MAID 2003 model 
and the dashed line is with the Lo+ multipole set to zero. The effect is particularly evident for Glt- Data are from flOllllLflgl and 
include statistical, systematic and model errors. 



TABLE 2. Background amplitudes 
for which the quadrupole amplitudes 
show sensitivity for MAID 2003 using 
the criteria listed in the text and the 
combined Bates and Mainz Q 2 = 0.126 
(GeV/c) 2 data. 

E1+ vs. E\^,M\-, F\,Fi 
L1+ vs. l[^ 2 ,Lo+,-F5 



seven terms shown in Table|2j when varied, lead to shifted central values or larger error bars for the resonant multipoles. 

1/2 3/2 

Some, like the E 1+ with M l+ are shifted but not outside the error bars and with no increase in the error bar size and 
so are not considered significant. 

Effort was made to search for a set of criteria using the correlation coefficient and the change in ^ 2 /d.o.f. that would 
identify all of the parameters which Fig.|7]identifies as significant. The criteria for significance were an increased error 
and/or a shift in central value. Either indicates a significant effect on the resonant multipole determination. In order 
to make the criteria robust, other models were put through the selection process as well. In addition to the MAID 
2003 model, Sato-Lee, SAID, and DMT were all used. The best identifier of significant parameters turned out to be 
the single test of \r\ > 0.7. In almost every case, this alone identified all the significant parameters. The % 2 sensitivity 
would identify some of the sensitive parameters but not others. 




FIGURE 6. Sensitivity to ¥ t and s and p in MAID 2003 using data at Q 2 = 0.126 (GeV/c) 2 . As in Fig. [5] the solid curve is the 
full MAID 2003 while the dashed line is now MAID 2003 with all the F,s turned off and the dotted line shows the results using 
only the resonant multipoles. This shows the significant effect the background has on some observables and that the F; amplitudes 
can also have effects of similar size. Data are from fill [Til [l2Tl and include statistical, systematic and model errors. 



EFFECT OF BACKGROUND ON RESONANT AMPLITUDES 

In order to see the effect the sensitive background amplitudes have on the extracted results, the resonant parameters 
resulting from each four parameter fit were plotted in Fig. [7] The horizontal bar indicates the position and error of the 
three parameter fit. The background amplitudes identified as significant do have an effect on the extracted multipoles 
relative to the three parameter fit. For each sensitive background parameter, the error increases and in most cases the 
central values shift. What is also interesting is that the F;s have a significant effect. This indicates that many small 
amplitudes can combine to have a large effect. 

To try to quantify the effect of the various background amplitudes on the resonant amplitudes, the RMS deviation 
of the various four parameter fits was taken for each model and identified as the intrinsic model error. For some fits, 
the RMS deviation was small and so the average of the four parameter fitting errors was used instead. In both cases, 
an estimate of the intrinsic error in each model was obtained. The results are shown in Fig.|8]and Table [3] along with 
the average and RMS deviation of the three parameter fits (model-to-model error). The figure and table indicate that 
the model-to-model variation is about the same size as the intrinsic model error (specifically for M\ + and the CMR 
the model-to-model error is larger while for the EMR it is somewhat smaller). However, the new error determination 
procedure is able to use one model alone instead of comparing it with other models. Each models's error can be 
assessed independently of the other models. 



While looking to improve the fits, an exhausive search was performed of all combinations of the three resonant 
parameters and any combination of the 13 remaining parameters. No significant improvement was found for either 
Q 2 = 0.060 or 0.126 (GeV/c) 2 . 

It is time, then, to look beyond fitting the multipoles. It is possible to modify internal model parameters (form 
factors, coupling constants) which affect many multipoles simultaneously but in different ways. This may allow the 
models to fit the data better. However, this fitting most likely needs to be performed by the model authors. 
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FIGURE 7. Effect of sensitive parameters on extracted resonance values using MAID 2003 for the combined Mainz and Bates 
data at Q 2 = 0.126 (GeV/c) 2 . Note the larger error bars and shifted central positions of certain combinations of multipoles. The 
horizontal band indicates the value and error resulting from the three resonant parameter fit. Similar plots were made for the 
remaining models. Their results are shown in Fig.[8]and Table[3] 



P 43 
§-42.5 
I 42 
-41.5 
f 41 

40.5 
40 

39.5 

r 0.5 

r "I 
1-1.5 
-2 
-2.5 
-3 
-3.5 
-4 

I -4 
c 

! -s 

-6 
-7 
-8 



i SL 

• MAID 
k DMT 



I | T SAID 

If i If { - RMSErr 





FMR 






i 








i ; tt 




• 

A 

■ 


* 


• - 



, CMR 




~ m 






E • » 


• T 


* 




i 


4 



FIGURE 8. Comparison of the estimated intrinsic and model-to-model errors for all four models at Q 2 = 0.126 (GeV/c) 2 , 
W — 1232 MeV. The intrinsic error was estimated using the RMS deviation across all the four parameter fits. If the deviation 
was small, then the average of the four parameter fit errors was used. Cases where the RMS deviation was used are indicated by 
the dashed vertical error lines. The horizontal band indicates the RMS deviation of the three resonant parameter fits using the four 
different models. 



TABLE 3. Preliminary fitting results and comparison of the estimated intrinsic and model-to-model errors 
for all four models for the combined Mainz and Bates data at Q 2 = 0.126 (GeV/c) 2 , W = 1232 MeV. See 
Fig.[8]for calculation details. 
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The understanding of the A will also be improved with experiments that are closer to complete. With target and 
recoil polarization, more observables are accessable and these have different combinations of multipoles. These new 
combinations will further constrain the models allowing better fits and smaller uncertainties in the backgrounds. Until 
new data are available, though, fitting the data and improving the models remain the best options. 

SUMMARY AND CONCLUSION 

Experimental results using the y*p — > pn° reaction have advanced the understanding of the shape of the proton and 
the A. However, the analysis process begins with extracting multipoles (which are not observables) from cross sections 
(which are). Without complete experiments including target and recoil polarization, the extraction must rely upon 
models for the background amplitudes. Performing standard three resonant parameter fits has allowed a good deal of 
progress to be made. Near resonance, fits using various models converge at Q 2 = 0.060 and 0.126 (GeV/c) 2 despite 
the differences in the model backgrounds. However, what has not been fully understood is the effect these differing 
backgrounds can have on the resonant parameters. 

To answer that question, we have added thirteen more background amplitudes to our three parameter fits and sys- 
tematically examined the effect of each one on all three resonant multipoles. Those additional background amplitudes 
are the four remaining s and p wave amplitudes, three isospin 1/2 amplitudes and six amplitudes we have constructed, 
the FjS. The large effect of some of the Fj terms shows how small multipoles which may have been ignored separately, 
can combine to have a sizable effect on the resonant amplitudes. 

As part of the systematic examination of the additional background amplitudes, correlations were found between 
them and the resonant amplitudes which led to larger errors and/or shifts in the values of the extracted resonant 
multipoles. We also found that while some amplitudes exhibit a large sensitivity in % 2 , no universal criteria could be 
found which would predict a sensitivity in the resonant amplitudes. Some background amplitudes which were sensitive 
did not affect the fits while others which were not sensitive did. 

However, varying the background amplitudes which were highly correlated with the resonant multipoles did affect 
the extracted resonant multipoles. Previous works have shown that the experimental and model-to-model errors are 
similar in size lfl8l \v% . What this exercise has shown is that the intrinsic model error is also similar in size to the 
model-to-model error. The current data really are challenging the existing models. So, without improvement in the 
models or more complete experiments, this is as far as the current data can take us. 

In general, the models agree with the data in a qualitative way but a good quantitative decription will require further 
refinement of those models. It is possible that some of the models may be made to fit the data much better with 
adjustment of the proper parameters. Adjusting a form factor or coupling constant within the model will change many 
multipoles in ways that are different from how they were varied in this study. Once the models are improved, they 
can be further tested with experiments that utilize target and recoil polarization. These introduce new combinations of 
multipoles which will further constrain the models. However, until new data are available, improvement of the models 
is the only option which will allow a better understanding of the A. 

Finally we return to the question that has primarily motivated this field: Do we have definitive evidence that the 
nucleon and A have non-spherical components and if so how large? Based on this study we follow reference J2l and 
present Figure|9]which indicates the final sensitivity to the quadrupole amplitudes. On the right is Gu which is sensitive 
to the Si_|_ quadrupole term. On the left is a special construction, Oe2 — °b(^o) + ^TriO^g) — °oifi%q = 7^121 1 101 1 1211 
which cancels out the dominant M\ + multipole contribution and shows the effect of the smaller E\ + quadrupole 
contribution. In Fig. |9l the range of predictions using all the four parameter fits was found by cycling through all 
the fits and storing the maximum and minimum. In this way, a high probability region was identified where the 
physical multipole would be expected to be. For comparison, the same procedure was repeated but with the quadrupole 




FIGURE 9. C£2 (left) and Glt (right) with error bands showing the range of all the MAID 2003 four parameter fits with (blue/dark 
region) and without (green/gray region) the quadrupole terms. The data are clearly in the region indicating quadrupole strength 
similar in size to that found in the models. These plots represent fits using the MAID 2003 model. The light region surrounding the 
MAID 2003 fit was found using the errors in the three parameter fit with fixed background terms. Data are from fToL fTTL fT^I and 
include statistical, systematic and model errors. 

amplitudes set to zero. On the basis of this study of the uncertainties in the resonance and background amplitudes we 
agree with the previous conclusion^ that a significant contribution of quadrupole amplitudes has been observed. 
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APPENDIX:CONTRIBUTION OF HIGHER PARTIAL WAVES IN THE LEADING 

MULTIPOLE APPROXIMATION 

The response functions can be expanded keeping only the terms which interfere with the dominant M\ + multipole. 
The multipoles for L>2 have been combined into the f)s in the following expansions which are called the Leading 
Multipole Approximation (LMA): 

4 MA = (^cos 2 e)|M 1+ | 2 

+ 6cos 2 0Re [Ef + M 1+ ] 

- sin 2 (3Re [E{ + M l+ ] + Re [M* l+ t\\ ) 

+ (1-3 cos 2 0) (Re [M* + Mi_] + Re [M* l+ Fb] ) 

+ 2cos0(Re[£ o * + M 1+ ] +Re[M* l+ F\]) (8) 



R^ A = -Re [M( + (3£ i+ + 37^+7^+ 3Mi_)] sin 2 
3 
2 



^|M 1+ | 2 sin 2 (9) 



R^ A = sin0Re[Mi + (7^ + Lo+ + 6Li + cos0)] (10) 

R L ^ = sin0Im[M 1 * + (7% + L o+ + 6Li + cos0)] (11) 

7 E2 A = -12sin 2 0Re[£ 1 * + Mi+] 
- 2Re [M* l+ 7^\ sin 2 

+ (2cos 0+2) (Re [£ *+ M i+] + R e [M^K] ) (12) 



